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We prove two statements about the long time dynamics of integrable Hamiltonian systems. In classical 
mechanics, we prove the microcanonical version of the Generalized Gibbs Ensemble (GGE) by mapping it 
to a known theorem and then extend it to the limit of infinite number of degrees of freedom. In quantum 
mechanics, we prove GGE for maximal Hamiltonians - a class of models stemming from a rigorous notion 
of quantum integrability understood as the existence of conserved charges with prescribed dependence on a 
system parameter, e.g. Hubbard U, anisotropy in the XXZ model etc. In analogy with classical integrability, the 
defining property of these models is that they have the maximum number of independent integrals. We contrast 
their dynamics induced by quenching the parameter to that of random matrix Hamiltonians. 


The past decade has witnessed an unprecedented experi¬ 
mental access to global coherent dynamics of many-body in¬ 
teracting systems. As a result, a new area that could be called 
“far from equilibrium many-body Hamiltonian dynamics”, 
“coherent many-body dynamics” or “quantum quenches” has 
emerged. A major part of research in this area has focused 
on testing the GGE|[Tl|2l in various integrable models. GGE 
refers to a density matrix or, in the case of classical mechan¬ 
ics, a phase space distribution function 

( 1 ) 

where are a (complete in some sense) set of integrals of 
motion for system Hamiltonian H and Z is a normalization 
constant. Suppose the system evolves with H starting from a 
non-stationary state. The statement of GGE is that the infinite 
time average of an observable O coincides with its ensemble 
average with the density matrix pH. 

Most authors test GGE in quantum models without clari¬ 
fying their notion of quantum integrability. The latter how¬ 
ever is a tricky concept with no generally accepted defini¬ 
tion, making the quantum GGE conjecture essentially unfalsi- 
fiable. The notion of classical integrability on the other hand 
is unambiguous It). Eor this and other reasons, it makes sense 
to first understand the status of GGE in classical mechanics. 
We will see that the microcanonical version of GGE - Gener¬ 
alized Microcanonical Ensemble - is exact for a general clas¬ 
sical integrable Hamiltonian. In a parallel line of inquiry, we 
will prove GGE for a class of models that emerge from a re¬ 
cently proposed complete notion of quantum integrability. 

Generalized Microcanonical Ensemble (GME) in classical 
mechanics is the following phase space distribution: 

n 

p{p,q) = V~^Y[6{Hk{p,q) - hk), ( 2 ) 

k^l 

where q = (qi,..., g„) and p = (pi,... ,p„) are the gener¬ 
alized coordinates and momenta. Suppose the system evolves 
with an integrable Hamiltonian H {p, q) starting from a point 
(POj 9o)- Let hk = Hk{po, qo) be the values of its integrals of 
motion for this initial condition. The statement of GME is that 


the time average of any dynamical variable 0{p,q) is equal 
to its phase space average with distribution 

lim ^ [ 0{t)dt= [o{p,q)p{p,q)dpdq, (3) 
T->oo T Jq J 

where O (t) = 0{p{t), q{t)). Eq. (j^ also holds for inte¬ 
grable classical spin Hamiltonians H{{sk}), in which case 
Pk = cos and Qk = where 9k and (jtk are the polar and 
azimuthal angles defining the spin direction. Eq. 0 is valid 
for any number of degrees of freedom n, so one can take the 
limit n —> 00 on both sides. Moreover, we will argue that the 
limits n —r oo and T —r oo commute (a tremendous simplifi¬ 
cation) as long as the frequency spectrum of 0{t) is free from 
a certain anomaly near the zero frequency. 

As a first step towards a similarly unambiguous statement in 
quantum mechanics, we also analyze GGE in the framework 
of a rigorous formulation of quantum integrability ||5l. Sim¬ 
plest models that arise in this approach are type-1 or maximal 
Hamiltonians - general N linearly independent commuting 
N X N Hermitian matrices of the form H{x) = T + xV, 
where a; is a real parameter. Type-1 matrices represent blocks 
of various exactly solvable many-body models (such as ID 
Hubbard and Gaudin magnets) for certain sets of quantum 
numbers (total spin projection etc. jSil and also describe e.g. 
a short range impurity in a metallic grain0. We prove GGE 
is exact for any N and explicitly determine Pk in Eq. Q. The 
GGE density matrix for quenches of the parameter x turns out 
to be non-thermal for type-1 Hamiltonians. In contrast, if we 
choose T and V randomly, the post-quench asymptotic state 
is thermal in W —oo limit. This emphasizes the importance 
of a well-defined notion of integrability as naively one could 
claim N integrals of motion (e.g. projectors onto the eigen¬ 
states) in the random matrix example too. We also relate the 
non-thermal behavior to localization. 

A characteristic feature of type-1 and classical integrable 
systems is that in both cases the number of independent in¬ 
tegrals is the maximum allowed by the definition. Their dy¬ 
namics are constrained by the integrals apart from linear in 
time phases (angles) that cancel out upon time-averaging or 
dephase in the thermodynamic limit. As the result, the inte¬ 
grals of motion fully determine infinite time averages. The 
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situation when the number of conservation laws is apprecia¬ 
bly less than the maximum is unclear and we will not consider 
it here. A common belief is that in quantum exactly solv¬ 
able systems this number scales as the logarithm of the size of 
the total Hilbert space, which agrees with most existing con¬ 
structions of conserved charges (integrals). Further analysis 
however reveals additional integrals |l6l[T0l, which can be cru¬ 
cial in identifying the proper ensemble ifTTll . It therefore seems 
likely that some version of GGE or GME does generally hold 
in quantum mechanics if integrability is properly defined and 
all independent integrals are taken into account. 

GGE apparently holds for relatively simple quantum mod¬ 
els, such as ID hard-core bosonsJTl and Luttinger liquids ifT^ . 
but fails e.g. in the XXZ and attractive Lieb-Liniger mod¬ 
els for a seemingly reasonable choice of integrals iFt. |fT3l - 
[la. More generally, ifTSlfTSll argue that GGE fails for models 
with bound states and ca, that it reproduces global observ¬ 
ables only in models mappable to noninteracting uncorrelated 
fermions. Given that Eq. ([^ is a rigorous theorem in classical 
mechanics, a natural question is: what are the reasons for such 
failure? Eirst, the equivalence between GME and GGE can 
break down already on the classical level. The standard argu¬ 
ment to go from the microcanonical to the canonical ensemble 
requires that energy be an extensive property and interactions, 
roughly speaking, short rangedlTS). It is straightforward to 
extend this argument to the generalized ensembles, but then 
each Hk must have these properties, which is not necessarily 
the case. Eor example, classical Gaudin magnets are well- 
defined integrable models with long ranged H^. 

Another set of problems arise from relying on an incom¬ 
plete definition of quantum integrability. Main issues here are 
how we understand the independence (nontriviality) and com¬ 
pleteness of a set of integrals. Classical integrability, for ex¬ 
ample, requires n functionally independent integrals. If we al¬ 
low all integrals in quantum GGE, one can simply choose the 
projectors onto the eigenstates of the Hamiltonian. With this 
choice of Hk Eq. Q is equivalent to the diagonal ensemble 
for any Hamiltonian and the statement of GGE becomes tau¬ 
tological. One might object that projectors are nonlocal (not 
short ranged). Eirst, this is not the case in models with local¬ 
ized eigenstates. More importantly, while locality might be 
necessary for going from the microcanonical to the canonical 
ensemble, it hardly is a legitimate requirement in the definition 
of quantum integrability. Indeed, there is no such condition in 
the classical case and, moreover, there are quantum Hamilto¬ 
nians, e.g. quantum Gaudin and BCS models ifTTlfTOl that are 
nonlocal, but otherwise bear all hallmarks of integrability. On 
the other hand, we do not expect GGE to hold for an incom¬ 
plete set of integrals, at least the theorem Q does not. 

Both Eq. 0 and quantum infinite time average are invari¬ 
ant with respect to the choice of any nondegenerate Hermitian 
operator within a given integrable family as the system Hamil¬ 
tonian. At this point we loosely define an integrable family as 
the set of operators that share the same integrals of motion 
Hk- Eor example, in the usual construction of the conserved 
charges for Lieb-Liniger, ID Hubbard, and XXZ models one 


of the Hk is the Hamiltonian, while the rest serve as its inte¬ 
grals. We can alternatively designate any other Hk or their lin¬ 
ear combination as the Hamiltonian without modifying Eq. 0 
and the time average. Interestingly, there is a combination, 

Ef/3 =/3-'^/3fciTfc, (4) 

k 

for which GGE coincides with the Gibbs ensemble. Eq. 0 is 
similarly independent of the choice of a nondegenerate Hamil¬ 
tonian (see below) within the classical integrable family. 

To see the above invariance, note that the time average is 
given by the diagonal ensemble 0, 

(0(f))j = \im ^ [ {O {t))dt = V (5) 

T^oo J ^ 


where Cm are the coefficients in the decomposition of the ini¬ 
tial state I in) into the eigenstates, which are shared by all Hk- 
Similarly, conditions (in|7Jfc|in) = Tr pHk that determine fik 
are the same. Let us also note that Eq. 0 is useful for a 
macroscopic system when the thermodynamic and T —oo 
limits commute. Arguments we make below about the order 
of n —oo and T —^ oo limits in Eq. 0 apply here as well. 

In what follows we first derive GGE for the maximal (type- 
1) Hamiltonians and then prove Eq. 0 . The approach of 13 ]- 
m to quantum integrability addresses many-body Hamiltoni¬ 
ans at the level of blocks (sectors) stripped of all space-time 
and internal space symmetries. Basic objects sae N x N Her¬ 
mitian matrices of the form H(x) = T + xV, where a: is a real 
parameter (Hubbard U, anisotropy in the XXZ model, mag¬ 
netic field in Gaudin magnets etc.) We say that H{x) is inte¬ 
grable if it has at least one commuting partner H {x) = T+xV 
other than a linear combination of H{x) and the identity. An 
integrable family is a set of commuting, linearly independent 
Hermitian matrices Hfix) = Ti + xVi, the most general 
Hamiltonian in the family being H{x) = diHfix). Thus, 
we define integrability simply as the existence of integrals lin¬ 
ear in the parameter ll20ll . 

There is a natural classification of integrable families by the 
number of independent commuting matrices they contain. We 
say that H{x) is type-M when this number is N — M + 1. 
The maximum number of linearly independent, linear in x in¬ 
tegrals H{x) can have is N. Then it is a type-1 or, equiva¬ 
lently, a maximal Hamiltonian. It turns out that any type-1 
Hamiltonian can be parametrized by 3A^ numbers di, Ei, 7 i as 
H{x) = N~^ Yhi diHfix), where 






£i Ed 


P ^3 = l*)0’l + li)(*l’Pi = b')0'l>and |j) are the shared eigen¬ 
states of Ti- Conversely, given arbitrary 2N real di, and N 
complex 7 i, Eq. 0 yields a type-1 Hamiltonian H{x). 
Normalized eigenstates |Am) of Hfix) read 




A^(A 

m 


= 


iTfcl" 


^ (A,n 


(7) 
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where Am is any of the N real roots of the equation 

Er^ = -- w 

Eigenvalues of general type-1 H{x) = dkHk{x) are 

k 

For example, dk = £k yields an interesting Hamiltonian that 
describes a short range impurity in a metallic grain||9l 

Hi^{x) =x\-j){j\+'^e^\i){i\, ( 10 ) 


with eigenvalues = Xm- Here I 7 ) = 7»l*)- 

The number of integrals for type-1 is maximal and equals 
the dimension of the Hilbert space N. We can show that GGE 
is exact for any N by matching the diagonal ensemble, i.e. 

(Am|e-^'=''''"''(")|Am)=^|Cmr dD 

With the help of Eq. (j^ this becomes 

=ln|cm|"+lnZ, (12) 

, ^k 


with an explicit solution 


lik 



(13) 


An overall shift I3k ^ Pk + c adds a constant c to the RHS 
of Eq. in view of Eq. (j^. The solution for j3k is unique 
apart from the shift, i.e. an arbitrary choice of Z. 

To see how type-1 differs from a general H{x) of the same 
form T -f xV (which has no nontrivial commuting partners 
in our definition), let us compare parameter quenches from 
Xi = 0 to X/ in type-1 Hamiltonians to those for random T 
and V. Let the initial state be an eigenstate |i) of T. 

The diagonal ensemble for type-1 according to Eq. 0 is 


2 . 

Nl{Xm-ed^- 


(14) 


A natural choice for Sk are eigenvalues of a random matrix 
from the Gaussian orthogonal ensemble (GOE) with mean 
level spacing 5 cx N~^ and we also set | 7 ip = N~^. Eq. 0 
implies £k-i < Xk < £k- It follows that in N —i' 00 
limit c>c for A^ in an infinitesimal vicinity of £i 

and c>c N~'^ otherwise. Note that this indicates localization 
of the eigenstate |Am) of H{x) in the space of eigenstates 
of Tiim. Eor simplicity, we assume that the matrix ele¬ 
ment (Am|0|Am) = fo{Xm) is a smooth function of Am 
(Am becomes a continuous real variable when N —> c»). 
Eq. § then yields {0{t))t = Consider next the 

GGE with Hp = Hi{x), i.e. /3fc = fidik- Its eigenvalues 


are = (Am — The ground state corresponds to 

m = i and is separated by a gap cx N from the excited states. 
The normalized density matrix is simply pmm = 5mi and 
Tr {pO) = foisi) = This GGE is non-thermal for 

any type-1 Hamiltonian other than Hi, e.g. for Hi^{x). 

Now suppose V is also random and uncorrelated with T. 
Eigenstates of H{x) = T + xV and T decorrelate at x = 
0(l) ll22ll . so that Cm are components of a random vector. Av¬ 
eraged over narrow energy windows (|cmP)£; = at large 
N, which corresponds to the infinite temperature Gibbs distri¬ 
bution. Note that already in type-1 T (or V) is arbitrary. A 
natural choice of T is a GOE random matrix. However, any 
choice of T severely constrains V to ensure the existence of 
commuting partners. It is precisely this correlation between 
V and T that also makes the density matrix non-thermal and 
eigenstates localized. Therefore, even though for a random 
H(x) one can take Hk in Eq. 0 to be the projectors onto its 
eigenstates, it is of no consequence because it does not intro¬ 
duce correlations between V and T and so constructed GGE 
is just the Gibbs distribution, i.e. Hp = H{x). From this 
point of view, the statement of GGE is not that Eq. 0 repro¬ 
duces {0{t))t, but that Hp is distinct from H. Note that the 
exponential form of p in Eq. 0 is unimportant at finite N. It 
could as well be a different function of Hp and we still would 
be able to match the diagonal ensemble. It however plays an 
important role in showing that Hp H in N ^ 00 limit, 
even though type-1 Hamiltonians are nonadditive. 

Now we switch gears to classical mechanics to prove 
Eq. 0. There are two necessary conditions: (i) the level set 
of Hk{p, q) = hk = const is compact and connected and (ii) 
the frequencies a; = (wi,... ,a;„) of quasiperiodic motion 
with H{p, q) are incommensurate (see below). The first one 
is a standard assumption in the Liouville-Arnold theorem to 
show that the dynamics is confined to invariant tori. It means 
that the motion is bounded and, roughly speaking, integrals 
are properly chosen. Consider e.g. a ID harmonic oscillator 
2H — p^ + The level set H{p,q) = const (ellipse) 

is connected, so E is a proper choice and Eq. 0 holds. If 
we instead take H = {H — hi){H — / 12 ) with 2 > 0 as 
our integral, the manifold H{p, q) = const is not always con¬ 
nected. For example, H = 0 corresponds to two oscillation 
amplitudes Ai 2 . The time average of e.g. q^ is either 7l^/2 
or A 2/2 depending on the initial conditions, while the phase 
space average with p = V~^6{H) is always (Af + Al|)/4. 

Let us go from {p, q) to action-angle variables (/, ip). The 
only information the reader needs to be able to follow the 
proof of Eq. 0 is that this is a canonical transformation, 
action variables depend only on Hk and vice versa, and the 
manifold Hk {p, q) = hk corresponds to a unique set of val¬ 
ues Q!m of action variables Im- Hamiltonian equations of 
motion in action-angle variables are pk — {H{I),‘Pk} = 
dH{I)/dik = uik and ik = {H{I),Ik} = 0. It follows that 
cp{t) = ipQ+tu). The motion is characterized by n frequencies 
ujk (or n periods). Incommensurability means there is no vec¬ 
tor m 0 with integer components such that m uj = 0. The 







4 


frequencies are incommensurate for most initial conditions in 
a nondegenerate system, 7 ^ 0||4l. 

Eq. Q in new variables reads 

71 

p{I) = {2^)-^\{5{h-au), (15) 

fc=i 


where we took into account that ipk varies from 0 to 2tt in 
determining the normalization constant. Since the Jacobian of 
a canonical transformation is 1, Eq. 0 becomes 

( 0 (f)), = 10(1, ^)pil)dld^ = J (16) 


where on the RHS we supressed the dependence of O on con¬ 
stants a = (ai,..., an)- Eq. (16 1 is a known theorem in 
classical mechanics called “the theorem on averages’ ’0. We 
prove it somewhat differently. 0{(p) is periodic in each ipk 
with period 27r. Expand it in multiple Eourier series 


(17) 

m 


where the summation is over all n-dimensional integer vectors 
m. The time-dependence of O along any trajectory is 

0{t) = 0{ip{t)) = (18) 


The finite time average of 0{t) is 


1 

T 


0{t)dt = Og + 


m^O 


^2'KiTmLJ — l 


2TriTm ■ u> 


^2-Kim-(-pQ 


(19) 


Since m • a; 7 ^ 0 for m 7 ^ 0, the last equation implies 

( 0 (f))t = lim f 0{t)dt = ag. 
T-)-oo 1 Jq 


Next we evaluate the RHS of Eq. < [T^ 

n p 2 tz 


o{^) 


dtp 

(27r)" 




fc=i ■ 


dpk 

2tt ' 


PirirrikVk 


— ^ ^ ^m^miQ - - - — ^ 0 - 


( 20 ) 


( 21 ) 


This completes the proof of theorem Q. 

Note that the phase space average is unconditionally equal 
to ag, while the time average equals ag only when the fre¬ 
quencies are incommensurate. Take, for example, a 2D 
anisotropic oscillator 2H = p\ -\- + ojiqi + u; 2 g| and 

O — (<?! + The phase space and infinite time time aver¬ 
ages are {Af -f ^ 3)72 and (Aj + A 2 + A 1 A 2 cos ^ 3 ) 72 , 
respectively, where a is the phase shift between the oscilla¬ 
tors. The averages do not agree when uji = uj 2 and, moreover, 
the time average depends on an initial condition other than the 
integrals (on a particular trajectory on the torus). 


As an application of the theorem, consider the semiclassi- 
cal Dicke model that describes coherent spontaneous emission 
from a large number of atoms (superradiance) interacting with 
a cavity electromagnetic mode ll^[25]| 


H = 2eSz + ujbb -f g{bS- + bS+), (22) 


where y/^b = {p — iojq) is the classical counterpart of the 
harmonic oscillator annihilation operator and 5” is a classical 
spin of length S. We are interested in the time averaged ampli¬ 
tude of the bosonic mode {\b{t)\)t- There are two degrees of 
freedom and two integrals. Hi = H and H 2 =bb + Sz, so H 
is integrable. In general | 6 (f)| is an elliptic function. A brute 
force evaluation of the RHS of Eq. Q yields a ratio of two 
complete elliptic integrals. The system has a stable equilib¬ 
rium at Sz = —S,b = 0 and an unstable one at = S,b = 0. 
Especially interesting are initial conditions in the vicinity of 
the unstable point. In this case | 6 (f)| first grows exponentially 
and then turns into a sequence of secant pulses (solitons). Let 
Sz = S and |&(0)| = r^. Eq. ([^ now gives 


(l&WD* 


21n(7rr^72fi)’ 


= 25'- 


(w — 2 e)^ 
4^2 


, (23) 


where is the maximum of |6(f)| ll23l . 

Einally, we turn to the analysis of n —> c» limit in Eq. 0 . 
Since the set of integer vectors m is countable, we can rewrite 
Eq. ( [TSl ) as a single sum 

0 (") (t) = a[,”^ + Y (24) 

k 


(n) 

where we separated the zero frequency term and Uq — ag 
in Eq. 0 assuming the frequencies uj^ are incommensurate. 
The superscript n indicates the number of degrees of freedom. 
Time averaging first and then sending n to infinity yields 

lim lim — / 0^^\t)dt = lim Og”^ = alT' (25) 
n—fCJO T—^00 T Jq n—foo 

Now consider the opposite order of limits. In n —c» limit 
the frequency spectrum generally consists of continuum and 
discrete parts. Nonzero discrete frequencies do not contribute 
to the time average, so we drop them for brevity, 

i-b 

0°°{t)= lim (26) 

n—^oo I 

J a 

Whether the limits T — 00 and n —>■ 00 commute depends on 
the behavior of the function F(H) = near H = 0 . 

If F(H) is integrable, the integral in Eq. (j2^ vanishes (de¬ 
phases) as f —?► 00 by the Riemann-Lebesgue lemma. Then 

lim(0(")(<))^ = (0“(f))^ = ao-, (27) 


i.e. the two limits commute. This is the most likely scenario 
in classical integrable many-body models. In fact, it is not 
obvious if there are any reasonable counterexamples. 
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The relevant quantity at large n is 0°° (t) and if the limits 
do not commute, the theorem ([^ looses its predictive power. 
This happens if F(fl) oc S(fl) at small fl. flk that vanish 
in n —oo limit must have an anomalously large weight or 
density i'(fl), so that a dramatic redistribution from finite to 
zero fl occurs as n —oo and a delta-function emerges from 
the continuum. Then the integral in Eq. (27 1 does not vanish at 
f —?► oo and 0°° (t) -ft off, while the n —>■ oo limit of the time 
average {0^^\t))t is still ag°. In linear analysis around the 
ground state, 12^, are the excitation energies (normal modes). 
F{^l) cx (5(f2) means we cannot write their contribution to 
0{t) in the usual way as an ordinary integral with the density 
of states. This signals a macroscopic degeneracy of the ground 
state in the thermodynamic limit. Similar analysis applies to 
the order of iV —c» and T ^ oo limits in Eq. Q. Note also 
that in the absence of isolated nonzero frequencies (assuming 
F{n) is well-behaved) 0°°{t) —?> const at large time. Then 
the strong version|0 of GME holds. 

Eor example, consider quenches of the detuning w in the 
many-body generalization of the semiclassical Dicke model 


H = '^2eksl+ujbb ++bs'^) . (28) 

k=l k=l 

This is also a classical integrable Hamiltonian 1261 . It 
describes an s-wave BCS-BEC condensate of atoms and 
molecules in the mean-field approximation. A relevant ob¬ 
servable is the superfluid order parameter A{t) = gb{f). Eor 
infinitesimal quenches we know A(f) as well as individual 
spins Si{f) at all times and finite n. It is straightforward to 
show that n ^ oo and T ^ oo limits commute. In particular. 


|A(f)| =Ao- 


k=l 


Cfc cos2Hfef, 


(29) 


where (e^ — -f Ag up to corrections of or¬ 

der n~^ and we suppressed the superscript n. Coefficients 
Cfc cx n~ In n —i' oo limit the second term in Eq. 
becomes an integral with a smooth F{£) = c(e)i^(e) and 


|A(f)| 


= lim„ 


, An as f 


oo, see 


for de¬ 


tails. Moreover, the strong version of GME applies. However, 
for stronger quenches there is a regime where |A(t)| asymp¬ 
totes to a periodic function, i.e. discrete nonzero frequencies 
are present in the continuum limit. Then the strong version no 
longer holds. The theorem ([^ still works and we expect the 
limits n ^ oo and T ^ oo to commute as before, though the 
analysis of this case is more difficult. 

In conclusion, we pursued two independent threads in this 
paper. Eor classical systems, we have seen that the GME 
given by Eq. ([^ holds for integrable Hamiltonian systems for 
a proper choice of integrals as long as the motion is bounded 
and the frequencies of the quasiperiodic motion are incom¬ 
mensurate. We argued that the limits of infinite averaging 
time and number of degrees of freedom typically commute. 


For quantum systems, we have seen that the GGE hods at any 
N for maximal N x N Hamiltonians that emerge from a rig¬ 
orous notion of quantum integrability and determined jdk in 
Eq.g. The characteristic feature of these systems is that they 
have the maximum possible number of independent integrals 
making them analogous to classical integrable systems. In 
contrast to random matrix example we studied, the GGE den¬ 
sity matrix for quantum quenches in maximal Hamiltonians is 
non-thermal. We discussed the potential reasons of failure of 
GGE and argued that one the main reasons is an incomplete 
set of integrals as a consequence of an ambiguous notion of 
quantum integrability. 
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